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We present search results based on next-to-leading order predictions for the pair production of 
color-adjoint leptons at the LHC. Quantum effects are sizable, dominated by pure QCD corrections, 
and sensitive to threshold effects. We illustrate the stabilization of scale dependences and confirm an 
excellent agreement between fixed-order and multi-jet predictions for representative distributions. 
Finally, we examine the trademark collider signatures of leptogluon pairs. Based on the CMS 
leptoquark search we derive a mass bound of 1.2-1.3 TeV for charged leptogluons, significantly 
improving the constraints available in the literature. 
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I. INTRODUCTION 

Leptogluons are color-adjoint fcrmions with non-zero lepton number. We can categorize them as strongly interacting 
partners of the Standard Model leptons. Leptogluons may arise, for instance, as a manifestation of compositeness pQ. 
Proposals along these lines are motivated by the flavor and high-scale gauge structure of the Standard Model which 
suggests a relation between quarks and leptons. This can be explained, for instance, if both types of particles share a 
common substructure and emerge as bound states of fundamental constituents. These building blocks will be bound 
together by a new confining force. The onset of a confinement scale A will characterize the typical masses of the new 
heavy composites. If the basic quark and lepton constituents turn out to be charged under the Standard Model gauge 
group the new composite quarks and leptons will inherit these gauge charges and will thus be pair-produced through 
their ordinary couplings to gauge bosons. Alternatively, they can be produced in association with their Standard 
Model counterparts, via higher dimensional operators which govern the radiative transitions between both. 

Historically, the seminal ideas on compositeness [3] cristallized into a paradigmatic class of models, the so-called 
preon models [ 4 j. With the term preon, as coined by Pati and Salam, one denotes these more fundamental building 
blocks of Nature, whose condensates make up the "elementary" quarks and leptons, together with a tower of excited 
states. For example, in the framework of the so-called fermion-scalar models, the Standard Model leptons emerge 
as bound states from the combination of a fermionic preon and a scalar anti-preon I = (FS); if these constituents 
are charged under SU(S) C , this naturally leads to leptonic composites of the sort 3 <£> 3 = 1 © 8: the color-singlet 
state can then be identified with an ordinary lepton, while the color-octet lepton corresponds to a leptogluon, e.g. 
the color-octet electron e 8 . This argument illustrates how novel fermionic states transforming under higher SU(3) C 
representations serve as telltale predictions of compositeness models. Scalar or vector leptoquarks are another example 
of such heavy colored resonances More generally, TeV-scale strongly interacting sectors are typically linked to 
some kind of compositeness [BHS]. Their collider imprints, for instance the production of exotic objects such as e.g. 
excited leptons [10], quark sextets and decouplets have been gathering attention in the recent years. 

Phenomenological studies of leptogluon production are available for a diversity of frameworks, ranging from ded- 
icated calculations of leptogluon pair production at e + e~, ep, and pp(p) colliders |12] to generic approaches where 
the characteristic signatures of exotic colored states are portrayed, not only at colliders [TTJ [13] but also for precision 
physics [13] . Very recently, studies on the prospects of leptogluon searches at the LHC [TS] and the LHeC [TS] have 
been published. 

Experimentally, the current 3er, model-independent mass bounds from collider searches yield m eg > 86 GeV for a 
stable charged lepton octet; and m v% > 110 GeV, in the case of its neutral analoge v% [17] • Dedicated analyses by 
JADE and HI more than two decades ago have extracted leptogluon mass bounds attached to specific compositeness 
scale choices, typically in the range of m e ± ~ 100 — 200 GeV for A ~ 1 — 2 TeV [18] . 

In this paper we analyse the production of leptogluon pairs to next-to-leading order (NLO) QCD. We compute the 
total leptogluon pair production rates and examine in detail the different quantum effects. The reduced theoretical 
uncertainties we illustrate via the stabilization of the NLO predictions with respect to variations of the (unphysical) 
renormalization and factorization scales. We compute representative kinematic distributions to NLO and compare 
them to independent results based on multi-jet merging. All these NLO calculations we carry out with the help of the 
fully automated package MadGolem 19 22 . Finally, we explore the characteristic collider signatures of leptogluon 
pair production and discuss the prospects for leptogluon searches at the LHC. 

II. LEPTOGLUON PAIRS TO NLO 

We analyse the phenomenology of leptogluons at the LHC assuming a minimal extension of the Standard Model. 
We entertain the existence of an additional generation of lepton octets which consists of one electrically charged, color- 
adjoint (Dirac) lepton — the electron /positron octet — and its neutral companion, a neutrino octet v^. Generically, 
we shall label these fields as Zg = e^,v%. For our analysis, no assumption needs to be made on their possible weak 
charges. For concreteness, and unless otherwise stated, we shall assume both states to be mass degenerate. This setup 
provides a minimal, but representative framework in which to explore the trademark collider imprints of leptogluons. 
Interestingly, neutrino octets from compositeness scenarios and gluino fields from supersymmetric models share the 
same quantum numbers (barring the lepton number). In this sense, neutrino octets will portray very similar aspects 
as gluinos, with no connection to the supersymmetric nature of the latter. 

Leptogluon interactions at the LHC are entirely driven by QCD. As SU(3) C adjoints leptogluons couple to gluons 
through the covariant derivative = (<9 m )ab + ig s {T c A^)ab, where denotes the gluon field, g s the 



Figure 1: Feynman diagrams for the pair production of leptogluons to leading order. We represent both partonic subprocesses: 
quark-antiquark annihilation qq — > and gluon fusion gg — > Igls- 

strong coupling constant, and T the adjoint SU(3) C generators. We write these generators in terms of the structure 
constants T A C = -if ABC . The Lagrangian that describes the gl&l& coupling then yields 

C D i> u (0 - m h ) ^ a D -g s f ABC ^ ^ ^ A° . (1) 

There are two partonic subprocesses which contribute to pair production of leptogluons at the LHC: quark-antiquark 
annihilation and gluon fusion (cf. Figure [lj : 

pp(qq) -> kk and pp(gg) -t l s k , (2) 

where l$l 8 can either stand for the production of a charged eg e$ or a neutral v s v s leptogluon pair. Throughout 
our analysis we set the central renormalization and factorization scales to the average final state mass /i° = /i_R,F = 
mi 8 , which leads to stable perturbative results [23]. [24] . Unless stated otherwise, we fix the LHC energy to \/S = 
8 TeV, while the default leptogluon mass choice is m; 8 = 1 TeV. As for the parton densities we use CTEQ6L1 and 
CTEQ6M |25) with five active flavors. The corresponding values for the strong coupling constant ot s are computed 
with consistent values of a s (/Xij) obtained via two-loop running from Aqcd to the required renormalization scale /jr, 
again in the five active flavor scheme. 

A. Production rates 

We first concentrate on the total cross section for leptogluon pair production, <j{pp — » computed to leading 

and next-to-leading order with the associated K factors. A comprehensive numerical survey we document in Table [T] 
We consider leptogluon masses ranging from 500 GeV to 1.5 TeV. From here on we shall concentrate on the production 
of electron-octet pairs, pp — » e^eg". The neutrino octet field vs only contributes as a new colored virtual particle, 
entailing a mass-suppressed and hence numerically mild contribution. The constant ratio 

gjpp^eje^) = ^ 
a(pp -> v 8 v&) 

reflects the fact that, so long as m e ± = m U8 , the differences between the charged and the neutral channels reduce to 
a symmetry factor 1/2, which accounts for the two identical particles in the final state. The results we present 
for three LHC nominal energies: = 7, 8, 14 TeV. 

The total cross sections are sizable and very strongly dependent on the leptogluon mass. Leptogluons with m/ 8 ~ 
500 GeV would be typically produced at a rate of 1 — 10 pb. TeV-scale leptogluons, in turn, range around O(10) fb. 
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Table I: Total cross sections o(pp — > efe s ) and corresponding K factors, for different leptogluon masses and LHC energies. 
The renormalization and factorization scales we fix to the central value = fip = fx — mi s . 
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The size of the QCD quantum effects we quantify as K = cr NLO /er LO . This may attain very large values, spanning 
the range of 1.8 to 2.0 for \/~S = 14 TeV, or even higher values 2.0-2.2 for \/~S — 8 TeV. These uncomfortably large 
if-factors we should not interpret as a breakdown of perturbation theory. Instead, they signal artificially low LO 
production rates for 0(1) TeV mass particles, a known problem of the CTEQ parton densities [21] [22]. The NLO 
cross sections are actually perturbatively stable. On the other hand, the growth of the K -factors as the leptogluon 
mass increases we can trace back in part to a non-trivial threshold behavior from the real and the virtual NLO 
contributions, which partially overcomes the phase space suppression of the NLO contributions. We will expand on 
all these aspects below. 

The total rates drop nearly three orders of magnitude when sweeping the leptogluon mass range from 0.5 to 1.5 TeV. 
A closer view into this mass dependence we display in Figure [2] Herewith we profile the cross section (top panel) and 
-K'-factor (down panel) as a function of m e ± for the different partonic subchannels. The LHC energy we fix to 8 TeV 
(left) and 14 TeV (right). The dominance of the gg-initiated mechanism for low-mass leptogluons — roughly one 
order of magnitude above qq — follows from the large gluon luminosity; for larger m/ 8 values this effect is inverted. 
The qq mechanism includes not only the contribution from the valence quarks it, d but also the suppressed flavor 
excitations from the second generation, which amount to less than 10% of the qq —> e$ e§ budget. The quark-gluon 
crossed channels qq, which arise purely at NLO, stagnate at the per-mil level given our treatment of the collinear 
divergences [22] . 

In more detail, we can read off Figure [5] that cr(gg — > IgTg) to LO depletes faster than a(qq — > Ish) with increasing 
leptogluon masses, i.e. the associated K factor K gg exhibits a steeper growth. Unsurprisingly, we encounter a similar 
behaviour for gluino pairs [26] . These features we can understand, on the one hand, from the scaling at threshold, 
as we will explain in detail later. Moreover, they can again be related to the respective parton luminosities. Heavy 
leptogluon masses probe larger values of the Bjorkcn x-variable - this being the region where the quark parton densities 
become more competitive as compared to the gluon ones. 

It is illustrative to compare the leptogluon production rates to similar QCD-driven heavy particle production, in 
particular squarks [2"2"H2"4"j . gluinos ^[25], and sgluons [21]. Based on MadGolem [H] and Prospino [251 124] we 
can establish the hierarchy 

a{pp Ish) ~ cr{pp -> gg) ~ O(10) X a{pp GG*) ~ 0(100) x a{pp ->■ qq*) . (4) 

The differences we may track down to larger color charges involved in the production of color-adjoints as compared 
to fundamentals and to the spin representation of the final-state particles. The production of fermion pairs like 
leptogluons and gluinos is more efficient than that of sgluons. Finally, it is also worth pointing out that the phe- 
nomenological profile of leptogluon and leptoquark pair production is very similar, in spite of their different spin and 



color representations. We will further exploit these coincidences in Section III 
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Figure 3: Sample Feynman diagrams describing the real emission corrections for the pair production of leptogluons via quark- 
antiquark annihilation (left) and gluon fusion (right). 

B. Real emission 

Real emission is part of the QCD quantum effects to leptogluon pair production. It comprises the O(a^) three- 
particle final states with a leptogluon pair plus an additional parton, which accounts for the QCD jet radiation, either 
from the initial-state partons or the final-state leptogluons (cf. Figure [3]). 

These contributions arc infrared divergent, as a result of the soft and collinear singularities developed by the initial- 
state radiation, as well as the soft poles from the gluon radiation off the leptogluon legs. These infrared divergences we 
subtract using massive Catani-Seymour dipoles [37J HI] ■ O n top of the Standard Model initial- initial and initial-final 
dipoles MadGolem automatically provides the corresponding leptogluon dipoles which are required for new final- 
final and final-initial singularities. Since the leptogluons are color octet massive fermions, their dipoles are equivalent 
to the massive fermions with the replacement of the color factors Cf — ► Ca- 

The MadGolem implementation is based on an extended version of MadDipole (2!Jj and retains a variable phase 
space coverage in terms of the FKS-style cutoff < a < 1 [30]. As the default value we use a = 1. The numerical 
performance of the leptogluon dipole we can assess from Figure [4] Notice that while the individual contributions 
from the integrated and the non-integrated dipoles diverge logarithmically for small a their sum is stable over roughly 
eight orders of magnitude. These results we compute for the uu — » e%e%g partonic subprocesses, assuming 1 TeV 
leptogluon final states. 



C. Virtual corrections 



The second class of NLO quantum effects includes 0{a 3 s ) terms from the exchange of virtual gluons, quarks and 
leptogluons, giving rise to self-energy, vertex and box one-loop corrections. A sample of the corresponding Feynman 
diagrams we display in Figure[5] Again, we compute them using MadGolem [IS], generating the one-loop amplitudes 
with Qgraf[3T]. To handle the helicity and color structures we rely on a dedicated set of currently in-house routines, 
prior to finally reducing their tensor structures in the Golem 32, 33J framework. The resulting one-loop integrals 
we numerically evaluate with the OneLoop library [33]. For a detailed account on the structure and performance of 
the tool, we refer the reader to Ref. [TU] . 

All divergences we regularize using standard dimensional regularization, analytically extending all integrals and 
internal propagators to n — 4 — 2e dimensions. We then subtract the infrared poles by including the integrated 
Catani-Seymour dipoles. The collinear higher-order corrections we absorbe into the parton densities. The ultraviolet 
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Figure 4: Dependence on the FKS-like a parameter corresponding to the final-final leptogluon Catani-Seymour dipole for the 
partonic subchannel uu — > e^e^g. The leptogluon mass we fix at m e ± = 1 TeV. 





Figure 5: Sample Feynman diagrams describing the NLO virtual corrections to the pair production of leptogluons via quark- 
antiquark annihilation (left) and gluon fusion (right). 

divergences in turn cancel against the corresponding ultraviolet counter terms, which are part of the leptogluon-nlo 
model implementation and are generated automatically. Its structure is completely defined by the on-shell renormal- 
ization of the leptogluon mass and the MS renormalization of the strong coupling constant, with a decoupling of the 
heavy colored degrees of freedom |35j . We check the cancellation of divergences both analytically and numerically. 
The finite parts for the one-loop amplitude we compare to an independent implementation of the leptogluon model 
in FeynArts, FormCalc and LoopTools [36 . In the Appendix we supply further details on the renormalization. 



We examine in detail the anatomy of the virtual corrections in Figures 6][7 The patterns share many common 
features with sgluon |21j and gluino |22 ( l26 j pair production. To start with, in Figure[6]we survey the mass dependence 
of the real and the virtual corrections for both parton channels gg and qq. For each of them we superimpose the total 
LO and NLO rates, together with the partial NLO virtual and real contributions <j virtual and cr rcal (upper panels), 
while the relative size of the NLO corrections we depict below. 

The different subsets of one-loop topologies we examine in Figure [7j Their relative weight with respect to the 



LO results we describe via the ratio Aa/a 



LO 



NLO.i 



a )/c , where i runs over all the different one- loop 



contributions — including gluon emission which enters through the integrated dipoles. Again, we separate the gg and 
qq subchannels. The crossed channel qg does not develop any virtual corrections and hence it is not included. Notice, 
though, that it is required to achieve a complete cancellation of the collinear divergences. The figure unveils, first of 
all, large vertex corrections to the QCD coupling ge%e% of at least 0(40%) and with opposite signs for each of the 
subchannels. The corresponding mass dependence turns out to be relatively mild, with a total variation not larger 
than 10% for the whole leptogluon mass range under scrutiny. 

Sizable (positive) corrections are driven by the integrated Catani-Seymour dipoles, i.e. from soft/collinear gluon 
emission. Also large are the box contributions for the gg initial state. Finally, leptogluon self-energy corrections 
become large for the gg initial state. Interestingly, they all feature a characteristic growth with the leptogluon mass. 
Both, the size and the acute mass dependence of these subsets we can trace back to genuine kincmatical and dynamical 
features. 
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Figure 6: Total LO and NLO leptogluon pair cross sections <r(pp — > e^e^T) as a function of the mass. We separate gluon fusion 
(left) and quark/antiquark annihilation (right). In the lower panels we account for the relative size of the corrections, separated 
using Catani-Seymour dipoles with a — 1 and including the integrated dipoles as part of the virtual contribution. 
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The steeper growth of the K factors for the gg-initiated subprocess can be in part understood from the LO threshold 
behavior: there, the rate scales like a gg ~ irmf^/S versus a q q ~ irmf^fi/S, with j3 2 — 1 — 4mf /S, which means that 
° :L ' {99 ~^ Ws) decreases faster than a ho (qq — >• Igls) with growing leptogluon mass. The dominant virtual corrections 
benefit from two threshold effects. An enhanced Sommerfeld rescattering [37] describes the exchange of soft gluons 
between slowly- moving final-state leptogluons, leading to a Coulomb singularity cr virtual ~ Tra s /(3. The complete 
logarithmic enhancement of the soft gluon emission behaves like a - Alog 2 (/3) + Blog(8/3 2 ). 

Both sources of threshold enhancements become more efficient close to threshold for the steep gluon luminosities 
and at large leptogluon masses, and can eventually be resummed [38] , 

The enhanced leptogluon self energies in the gg — > Z 8 ^s case can be correlated to the leptogluon mass insertion in 
the t and u-channel leptogluon exchange - with no counterpart for qq — > Igls- The remaining one- loop contributions, 
consisting of QCD-like quark-gluon and triple-gluon vertex corrections and the box and self-energy corrections for 
qq —> Ish are essentially featureless and contribute with overall yields below 10%. 



D. Scale dependence 

One of the main motivations of higher-order predictions is to stabilize the (unphysical) dependence on the renor- 
malization and factorization scales. We introduce them when we remove the collinear and ultraviolet divergences at 
a given order in the perturbative series. In the limit where all terms in such series are retained the scale dependence 




Figure 8: Leptogluon pair cross section cr(pp — > eg~e g ) as a function of the leptogluon mass. The scale uncertainty envelope 
corresponds to an independent renormalization and factorization scale variation in the range A t °/2 < (ir,f < 2/i°. 
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Figure 9: Renormalization and factorization scale dependence for leptogluon pair production pp — > eg eg . The plots trace a 
contour in the /ir-/if plane in the range fj, — (0.1 — 10) x fj? with fi° = m ±, and m e ± = 1 TeV. 



vanishes. In practice, we truncate the expansion on the coupling constant at a certain order C(a"). In the absence 
of any systematic correlation between different orders n we expect the scale dependence to englobe the asymptotic 
cross-section values. This expectation has been well confirmed for many QCD-mediated processes [2TTf24] , 

In Figure [8] we depict the evolution of the total leptogluon pair cross sections at LO and NLO with the leptogluon 
mass, for the LHC operating at 8 and 14 TeV. The band envelope accounts for the theoretical uncertainty, which 
we estimate as conventionally by varying the scale choices within the range fi /2 < hr.f < 2/x°. We see that over 
the entire mass range the LO and NLO curves hardly overlap, which means that at least for the LO cross section a 
scale variation by a factor two around the central scale would not have been a conservative estimate of the theory 
uncertainty due to missing higher orders in perturbation theory. 

Complementary vistas on the stabilization of the scale dependence we provide in Figure [9] We delineate the LO 
and NLO cross sections for an independent variation of the factorization and renormalization scales. The five panels 
correspond to a contour in the {1r — fip plane, which we define in the left panel. The smoothing of cr(fi) at NLO 
nicely illustrates the stabilization of the higher order prediction with respect to scale choices. Quantitatively, the LO 
variation stays within Aa LO /a LO ~ 0(65%), while at NLO it reduces to Acr NLO /cr NLO ~ 0(30%). 



E. Distributions vs jet merging 



Relying on fixed order NLO predictions we can extract a suitable normalization to the event rates from standard 
Monte Carlo simulations. Nevertheless, one still needs to assure that this framework can be safely promoted to 
the main distributions. It has been shown in the literature that transverse momentum and rapidity distributions 
are relatively stable with respect to higher order corrections for heavy particle production [20-23]. As long as the 
collincar approximation includes sizable p? values owing to even larger final-state masses QCD jet radiation should 
be properly accounted for by the parton shower description [3HI 20] ■ 



In Figure 10 we check these expectations. Therewith, we display the transverse momentum and rapidity distributions 
for 1 TeV leptogluons. The different LO and NLO pieces we display separately. Real emission and virtual NLO effects 
we disentangle by introducing Catani-Seymour dipoles and separating the corresponding phase space regions within 
a = 10~ 3 . We can then compare the fixed-order NLO results from MadGolem with a matched matrix element and 
parton shower description using the MLM scheme including up to two hard jets [41] [42] . For the latter we employ 
MadGraph5 [101113] interfaced with Pythia [33]. We explicitly check that including just one hard jet barely changes 
the results, in line with the general statement that in our case the collinear parton shower is not a bad approximation. 

The two main histograms we normalize to unity, with the different NLO contributions shown to scale. Both 
approaches display very similar shapes, featuring central leptogluons with transverse momentum peaking at Pt(ss ) — 
300 GeV. Only when analysed individually, the px curves for the real and virtual NLO corrections depart significantly 
from the combined curve. The observed shift towards slightly harder and more central leptogluons in the merged 
result, which also appears for sgluon pairs [3T], we can deem to the additional recoil jets. 

Finally, in Figure 11 we present a panoply of Figure 10 now displaying the transverse momentum and rapidity 
dependence of the iVfactors including the scale dependence envelope around the central scale fj, = m ± . Remarkably, 
the if -factors remain relatively constant within the remaining NLO scale dependence as long as we stay within central 
rapidities and pt values around m e ±/2 > 500 GeV. This justifies the use of the if-factor as an overall reweighting 
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Figure 10: NLO distributions for leptogluon pair production as a function of the transverse momentum pr(eg ) and rapidity 
y(e% ). For the NLO curves we separately display the LO, virtual, and real contributions, separated using Catani-Seymour 
dipoles with a — 10 -3 . In addition, we show the distributions based on LO matched samples with up to two hard jets. The 
NLO and matched curves we normalize to unity while the different contributions to the NLO rates we show to scale. 
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Figure 11: K factor as a function of the leptogluon transverse momentum pT(e^) and rapidity y(e^). We include the scale 
uncertainty envelope for an independent renormalization and factorization scale variation within /i°/2 < (j,r,f < 2fj,°. 

factor for the leading-order event samples as long as we do not probe extreme phase space regions - see e.g. [35] for 
a counter example. 

III. LEPTOGLUON SIGNATURES 

From here on we will focus on trademark collider signatures of leptogluon pairs and discuss them in the light of 
the current LHC searches. Specifically, we discuss how to make use of ongoing CMS searches for leptoquarks [3PJ to 
provide updated leptogluon mass bounds. 

Besides the SU(3) C gauge couplings leptogluons can also interact with matter fields. These interactions are ulti- 
mately induced by the dynamics of the underlying model. They appear as non-renormalizable couplings which we 
can parameterize through higher dimensional effective operators. For simplicity, we limit ourselves to the dimension-5 
chromomagnetic operator. It is a representative example of such non-renormalizable couplings and accounts for the 
radiative decay of a leptogluon decay into a lepton and a gluon. This process is governed by the effective Lagrangian 

£ = TK^l ^ G U a ^ p t + a * p R) ^ + h - c - . ( 5 ) 



10 




Figure 12: Sample Feynman diagrams describing the characteristic dijet-dilepton final-state from leptogluon pair production 
followed by a radiative transition into a Standard Model lepton and a gluon each. 



with A denoting the characteristic scale of the ultraviolet completion, e.g. the compositeness scale of a hypothetical 
preonic model. The spinors tpf^ and tpi describe the leptogluon and Standard Model lepton. In the presence of the 
gluon we factorize the strong coupling constant g s . The overall coupling strength is affected by the parameters cll/r 
which we expect to be of 0(1). If instead of a strongly-interacting ultraviolet completion this transition is mediated 
by quantum effects we expect another suppression factor 1/167T 2 . Henceforth, we assume cir = 0, in agreement with 
electron chirality conservation and the left-handed nature of the Standard Model neutrinos. 

In the decay of leptogluons into a Standard Model lepton and a gluon both (massless) daughter particles emerge as 
back-to-back hard objects in the leptogluon rest frame. The leptogluon partial width we compute from the effective 
Lagrangian Eq.([5]), and we may cast it into the following form 

„ 9s m l 8 , 2 , 2 \ (n\ 

167r# (i+ r) ' () 

This form coincides with the width of a gluino decaying into a gravitino-gluon final state |47j after matching the 
respective effective couplings. Such analogies between leptogluons and gluinos are also rubber-stamped into their 
respective collider footprints, cf. Ref [35J for a detailed account on associated gluino-gravitino production at the LHC. 



We show a sample of leading-order Feynman diagrams for the pair production of leptogluons and their subsequent 
radiative decay in Figure [12] The overall rates we can describe within the narrow width approximation, 

a{pp -> l&h -> llgg) ^ a(pp -> kh) x [BR(/ 8 -> lg)} 2 . (7) 

The narrow width approximation is manifest from Eq.([6]) and assuming 4(mA 2 « 1. As a consequence the lep- 
togluon signatures will essentially not depend on the effective higher-dimensional interactions modelling the leptogluon 
decay. We should be able to derive universal constraints on the leptogluon masses unrelated to the characteristic ul- 
traviolet scale A and the coupling strength auR) of the effective l%lg interactions. 

As already mentioned, searches for (charged) leptogluons at the LHC should close follow searches for leptoquarks, 
as both processes share the same final-state signature [2B]. It consists of two isolated hard leptons alongside two 
central jets. These leptons get replaced by neutrinos or missing energy when we instead entertain the case of neutral 
leptogluons. As heavy resonances, leptogluon decays naturally give rise to hard px profiles for the resulting products, 
as well as sizable invariant masses for the dilepton m(e + e _ ) and leading-p^ dijet systems m(jj). Irreducible Standard 
Model backgrounds are Z/j* /W+jets, W + W~ +jets, and tt production. Additional subleading sources are single top 
production, electroweak dibosons, and QCD multijets, typically well below the per-cent level. 



A. Signal vs background 



The signal and backgrounds in this section we simulate with MadGraph5 [33] • A minimal extension of the 
Standard Model including the leptogluon fields and corresponding interactions we implement via Feyn Rules [49], 
which also provides the corresponding UFO model files [50]. We use MadAnalysis5 [5T] to analyse the event 
samples. The effective coupling/scale ratio a^/A we fix to 10 -5 , with clr = 0. While this choice has barely any 
impact on the total 0(a s ) leptogluon rates it suppresses additional 0{a s a 2 L ^ R / K 2 ) channels from gluon fusion and a 
i-channel leptogluon exchange. For the leptogluon mass the default value is mi 8 = 1 TeV. For the LHC energy we 
assume \/S — 7 TeV. The leptogluon effective couplings are diagonal in generation space, such that we search for 
electrons and electron-neutrinos, I = e^ 1 ,v e . 
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Figure 13: Correlated distributions for charged leptogluon pair production giving rise to 2j + 2l final-states (top panels) and the 
leading Standard Model background pp — > Z+ jets (bottom panels). We explore the correlations between the hardest lepton 
and jet transverse momenta (left) and the invariant masses of the dilepton and the hardest jet-lepton systems (right). 



As this process involves the production of colored heavy fields, we expect a sizable influence from the initial and 
final-state QCD radiation. This gives rise to multi-jet configurations which depart from the mere LO expectations and 
leave us with the separation of jets into QCD radiation and decay jets [HJ 301 H3 [52] . We model QCD jet radiation 
using MLM jet merging as implemented in MadGraph5 to generate signal samples including up to two extra hard 
jets. Technically, we enforce m™ m = wi™ m = 10 GeV, to avoid phase space singularities in the backgrounds. 

Representative signal and background correlations between different distributions involving the hardest lepton and 
jet decay products we examine in Figure |13[ They nicely illustrate the relevant kinematical features which we can 
rely on for a suitable search strategy. A first telltale imprint of the charged leptogluons involves hard charged leptons, 
each of them stemming from the decay of one of the leptogluons. This translates into very large dilepton m(e + e~) 
and dijet m(jij2) invariant masses. The natural scale for the signal selection is m(e + e~) > rriz- The leptogluon mass 



reconstruction, in turn, should proceed through the lepton-jet invariant mass m(ej) which we also display in Figure 13 
Another hallmark involves the transverse momenta of the hardest jet and the outgoing charged leptons. Additional 
jet emission from the parton shower is responsible for a smearing of this peak. Since this feature is correlated with 
the leptogluon mass it will depart very visibly from the background structure, as shown in the bottom-left panel of 
Figure [13] 

In analogy to the charged leptogluon case we show the relevant correlated distributions for the neutral leptogluon 
signal pp — > v$ v% — > ^ T +jcts in Figure 14 In the absence of identified leptons the experimental prospects are more 
challenging. After using a lepton veto to suppress the ti and (Z/j* — > e + e~)+jets backgrounds the leading irreducible 
background is Z-\- jets production with an invisible decay Z — ¥ vv. The key observable is missing transverse energy 
The expected $t from Z+jcts production displays a much softer profile than the signal. We illustrate these 
features in Figure [14] 

More generally, signal and backgrounds substantially deviate from each other in the number of final-state jets. 
Owing to the SU(3) C adjoint charge of the leptogluons as well as to their large masses we expect the leptogluon events 
to yield larger jet multiplicities, as reflected in Figure [15] By the same token, we also retrieve a much harder total 
visible energy Ht = Y) j +pt(s + ) +_Pr(e~), whose peak follows once more the leptogluon mass. 
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Figure 14: Correlated distributions for neutral leptogluon pair production giving rise to 2j + f T final-state signatures (top 
panels) and the leading Standard Model background pp — > Z+jets (lower panels). We explore the correlations between the 
leading and subleading jet transverse momenta (left) and the missing transverse energy vs the invariant mass of the leading 
dijet system (right). 



B. Data 



After establishing the NLO description of leptogluon pair production and studying the main kinematic features of 
the signal and background we are now ready to compare the precision predictions to LHC data. Specifically, we rely 
on the leptoquark searches by CMS |46j . Leptoquarks (lq) are described by a minimal extension of the Standard 
Model including one additional color triplet scalar. Our discussion henceforth will cover charged leptogluons only. 
The event samples for the signal are originally generated using Pythia 6.422 for a range of leptoquark masses 
rni q = 250 — 900 GeV, with a central rcnormalization and factorization scale /i — mi q . The major backgrounds are 
determined either from control samples in data or from Monte Carlo simulations normalized to data in selected control 
regions. The available experimental data corresponds to the = 7 TeV run with 5 fb _1 of integrated luminosity. 




Figure 15: Differential cross sections as a function of the jet multiplicity and the total visible (resp. missing) transverse energy 
for charged (left panels) and neutral (right panels) leptogluons. We also include the leading background pp — > Z+ jets. All the 
histograms we normalize to unity. 
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Figure 16: Distributions for the lepton-jet invariant mass (left) and the scalar transverse energy (right) for leptogluon pair 
production followed a decay into a lepton and a jet. The backgrounds we display as cumulative histograms and are taken from 
[461 153] . The leptogluon signal is reweighted to the NLO prediction. For comparison, we also include the leptoquark signal for 
the same mass. The bottom panels illustrate the relative differences between the leptogluon and leptoquark distributions. 



For further details we refer the reader to Ref. [46]. 

We now replace the leptoquark signal in the CMS leptoquark search with simulated leptogluon events, reweighted 
to the NLO rate for each assumed leptogluon mass. Following Ref. |2B], we model the triggers by a set of initial 
acceptance cuts. They are tailored to select potential e + e~ jj candidates among those events with two hard central 
leptons and QCD jets. Moreover, the jets are required to have a spatial separation of with respect to the electron 
candidates, while the dilepton invariant mass is enforced to fulfill 

p T (e) > 40 GeV |ry e | < 2.5 m(e+e") > 60 GeV Sf > 250 GeV 

p T (j)>30GeV M < 2.4 AR ej > 0.3 . (8) 

The cut on the scalar transverse energy — PT{e + )+PT{ e ~)+PT(ji)+PT(j2) significantly reduces the combinatorial 
backgrounds. 

After this selection the data is well compliant with the Standard Model background expectations, as reflected in 



Figures 16 and 17 Figure [16] shows the electron-jet invariant mass m e j (left) and the scalar transverse energy S*j> 
(right) after acceptance cuts. Data we display as black circles while the backgrounds form cumulative histograms and 
are taken from (3S] [S3] . Finally, we overlay the leptogluon and leptoquark signals (solid lines) , normalized to the 
respective NLO rates [M] . The masses of both resonances we fix to m = 400 GeV, to provide a replica of Figure 2 
from Ref. |46j including the leptogluon pair signal. In the lower panels we show the relative bin-by-bin differences 
between the leptoquark and the leptogluon predictions for identical masses. To deal with combinatorial issues when 
reconstructing the leptogluon mass, we distinguish two kinematic regions <j) G [0, ±7r] and we pair each of the leptons 
within a given event to the highest pt jet of the same hemisphere. 

Adopting the leptoquark analysis we can derive an approximate leptogluon mass. Our starting point is the 
acceptance-reconstruction efficiency factor as a function of the leptoquark mass, provided by CMS 53J. This we 
combine with the leptoquark acceptances to get the corresponding mass-dependent efficiencies e\ q . The acceptances 
we estimate resorting to a leptoquark simulation within MadGraph5 and implementing the cuts in Eq.([8]). The same 
strategy we follow to evaluate the corresponding leptogluon acceptances. Finally, we assume that the leptoquark ef- 
ficiencies eiq can be exported to leptogluons, ei q ~ e es , and derive the combined acceptance-reconstruction efficiency 
factors for leptogluons as a function of m eg . 

We observe that the discrepancies between the leptoquark and leptogluon acceptances, and so on the resulting 
acceptance-efficiency factors, are indeed small. In spite of the intrinsically distinct kinematics of color-triplet scalars 
vs color-octet fermions, non-negligible differences arise primarily in phase space regions hardly passing the cuts of 
Eq.Q. They mostly attain tails in the transverse momentum and rapidity distributions. We study these differences 
for a number of representative distributions and find them typically at the 0(10)% level. This is also indicative that 
eiq ~ e es is a reasonable hypothesis. 
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Figure 17: Expected and observed 95% CL limits on leptogluon pair production, followed by a decay into a lepton and a jet, 
based on the CMS leptoquark searches |46l 153] . The uncertainty bands represent 68% and 95% CL. We overlay the predicted 
production rates to NLO, including the theoretical uncertainty estimated from a full scale variation. 



Finally, we take the CMS leptoquark limits from a modified frequentist confidence level (CL) analysis [55] and 
translate them into the leptogluon limit, extending them to leptogluon masses m eg = 900 GeV — 1.5 TeV. This 
procedure is valid so long as the reconstruction efficiency does not change significantly within this mass range. This 
assumption is justified given the rather general analysis strategy and the experimental setup [53] , 



In Figure 17 we show the expected and observed 95% CL upper limits on the rate for charged leptogluon pair 
production as a function of the mass. The uncertainty bands around the median expected limits correspond to 1 and 
2a. We overlay a two-sided magenta band accounting for the theoretical uncertainties using the NLO predictions, 
derived from the independent variation of the renormalization and factorization scales around the central value, 
[i°/2 < fi < 2p°. The intersection of the expected and the predicted rates defines a 2a exclusion of m es < 1.2 — 1.3 TeV. 
This result constitutes a major improvement of the leptogluon mass bounds quoted at present |17j . 



IV. SUMMARY 



We have studied the pair production of color-octet leptons at the LHC. These leptogluons constitute a minimal 
extension of the Standard Model by charged e§ and neutral isg fields coupling to quarks and leptons via higher- 
dimensional operators. 

In the first part of the paper we have reported on the complete NLO calculation of the LHC production rate, 
using MadGolem. The genuine new physics structures like Catani-Seymour dipoles and field as well as mass 
renormalization are part of the automated setup. We find large production rates around O(10) fb for = 8 TeV 
and leptogluon masses in the TeV range. Significant quantum corrections arise mostly from pure QCD, i.e. virtual 
gluons and the massless initial/final state radiation. Leptogluon-mediated effects are mass-suppressed. The remaining 
theoretical uncertainty reduces from ~ 65% at LO down to ~ 30% at NLO. The NLO predictions are stable across 
the relevant phase space regions and agree well with results using LO multi-jet merging. This allows us to re-weight 
event simulations using jet merging to the NLO total cross section. 

In the second part we have studied LHC signatures of leptogluon pair production. We have described the leptogluon 
couplings to matter fields in terms of a dimension 5, generation-diagonal l$lg chromomagnetic interaction. Leptogluon 
pairs appear as characteristic dijet plus dilepton signatures, featuring hard transverse momenta and large invariant 
masses. These quantities can easily be distinguished from the irreducible background, mainly governed by pp — > 
Z+jets. For charged leptogluons, the signature is the same to leptoquarks. 

Finally, we use this fact to promote the current CMS leptoquark searches to leptogluon searches. We exclude the 
existence of charged leptogluons in the above framework with masses below ~ 1.2 — 1.3 TeV. 

Aside from its phenomenological relevance our work shows how the automated MadGolem framework can be used 
to efficiently study new physics models at the LHC. It completely automatizes the calculation of NLO cross sections 



15 



and distributions for pair production of new particles. MadGolem is an independent modular add-on to MadGraph 
and can easily be interfaced with its multiple user options and analysis tools. Following the present final testing phase 
we will make the code available to the theoretical and experimental communities at the LHC. 
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Appendix: renormalization 



Basic setup The ultraviolet renormalization counter terms we generate automatically from the Qgraf leading- 
order amplitude. The respective field, mass and g s renormalization constants we conventionally phrase in terms of 
two-point functions. They are part of the leptogluon model implementation and supplied as a separate library. At 
present, MadGolem fully supports the calculation of NLO QCD corrections for the Standard Model, the MSSM and 
several other extensions of the Standard Model featuring heavy colored resonances (e.g. sgluons). All the necessary 
renormalization constants we define through the relation between the bare and the renormalized fields, masses and 
couplings - the former being denoted with a "0" superscript: 

*(°)^4 /2 * ^my+Smy gi°^g s + Sg s , (9) 

where \& stands for all the strongly interacting degrees of freedom of the model, in our case ^ = q,g, lg. 

Given a generic QCD interaction with Lagrangian density C(^,m^,,g s ), we define its associated counterterm 
<5 jC (\& , m$, g s , <5\l/, 5m<f, 5g s ) as 

£(* (0) , m ( °\gW) = £(*, m*,g s ) +SC{^>, to*, 3s , 5*, Sm*,6g s ). (10) 
For further details on our notation setup, we refer the reader to Appendix D of Ref. [22] . 



Renormalization of the strong coupling constant In the presence of the leptogluon field, the /3-function of the 
Standard Model, and thereby the strong coupling constant renormalization, get modified. Using standard notation, 
we may formally express the quantum corrections to the quark-gluon (qqg) vertex in terms of the gluon (Z3) and 
quark (Z2) field-strength renormalization constants: 

z 1 = z g z 2 zl /2 . (11) 

These factors we can expand to order 0(a s ) as Zi = 1 + 5, + O(a^), with the counterterms Si being conventionally 
defined in the MS scheme. The strong coupling constant renormalization at one loop we can thus write as 

5g s = S 1 -S 2 -^S 3 . (12) 

Lcptogluons only couple to matter fields through higher-dimensional operators, such as the leptogluon-lepton-gluon 
interaction described by Eq.([5|. Consequently, quark self-energies - and hence 82 - are not affected by the presence 
of the leptogluons. The fermionic contribution to the (3 function remains therefore unaltered with respect to that of 
the Standard Model: 

SF = ~C F A € . (13) 

The shifted pole in the MS prescription is A = = — je + log(47r). The same reason explains the absence of leptogluon- 
mediated corrections to the quark-gluon vertex at one-loop, which means that 

5 1 =Sf M = -^(C A + C F )A e . (14) 

47T 
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Figure 18: Feynman diagrams for the leptogluon field renormalization at 0(a s ) and for the leptogluon-mediated corrections to 
the gluon field renormalization at 0(a s ). 

Instead, the pure QCD gluon/leptogluon interaction furnishes a novel piece to the gluon self-energies, as displayed 
in Figure [18} The latter results into 



5 l i=-sC A 



3tt 



(15) 



with s = 1(1/2) for ig = ef(vg,). Adding up the above result to the Standard Model contribution we get 



3?r 



(16) 



The final expression for 5 g s we obtain by plugging Eqs.( 14|-( 16 1 back into Eq.(12|, while explicitly decoupling 



the heavy (H) colored states - in our case, these are the top-quark and the leptogluon field (s). We implement 
such subtraction in a conventional zero- momentum scheme [23] . as described in Ref. [35] . That way we leave the 
renormalization group running of a s to be merely determined by the light (L) degrees of freedom. This corresponds 
to the definition of the measured value of the strong coupling, for example in a combined fit with the parton densities. 
The renormalization constant finally reads 



8 9 s = 



'4tt 2 



Air 



1 (m 

3 log b 



o C A 1 m 



(17) 



Both light (L) and heavy (H) colored particles contribute to the coefficient of the beta function, which in our leptogluon 
model reads fto = /3q + /9q , with 



0o 
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Ca — Cf N F T R 



and 



C F T R + As 



C A 



(18) 



By default, MadGolem sets the number of active flavors to Np = 5. Standard conventions for the SU(3) C color 
factors Cf — 4/3, Ca = 3 and Tr = | are assumed throughout. 



Renormalization of the leptogluon sector As a massive, color-adjoint particle, the leptogluon 2-point correlation 
function receives 0(a s ) corrections due to virtual gluon interchange, as shown in Figure 18 The corresponding UV 



divergences we can absorb into a redefinition of the leptogluon mass mi 8 and field-strength S Zi % . Both quantities we 
renormalize in the on-shell scheme. The renormalization conditions require i) the renormalized leptogluon mass to be 
the pole of the real part of the corresponding propagator; and ii) the renormalized propagator to have unit residue. 
These conditions we can cast in the following guise: 



lim 

p 2 —tmf 



ReE l8 (p)Vi 8 (p) = 0. 



(19) 



On the other hand, the (real part of the) renormalized leptogluon self-energy, 3?eS; 8 , we conventionally write as: 



3?e 



»Pr,Ez 



R t-'lsiRK 



TO/„E 



(20) 



where the scalar components L, R and S render: 



^h,L/R{ 



SZ, 



h,L/R 



2 m h 



(21) 
(22) 
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-g. f ABC < r [S9s + ± (8 z l8 + s z h + s z g )] A° 



iip, (f~mi 8 )SZi 8 ipi 8 - 8m l8 i/> h 



Table II: Counter term Feynman rules for the leptogluon-mediated interactions. 



Enforcing the on-shell renormalization conditions Eq.(19) on the above leptogluon 2-point functions, we can finally 



anchor the field and mass renormalization constants as follows: 



8Zi s ,l/r = -$te£ htL/R (mf s )-mt s Re [S' ig)L (mf g ) + S' /8)H (mf 8 ) +2S / i8 , s (mf g )] 



Smi 8 



-Ke [E, 8iL (t 



(23) 



Notice that the above expressions Eq.(20)-Eq.(23l equally hold for both charged and neutral leptogluon states. More- 



over, as far as we disregard electroweak effects, in practice we will find no distinction between the respective left and 
right-handed components. The conventional shorthand notation f'(p 2 ) = df(p 2 )/dp 2 is employed therewith. 

The analytic form of all renormalization constants we reduce down to 1 and 2-point scalar one-loop functions, in the 
conventions of |34j . UV divergences we handle by means of a standard 'tHooft-Veltman dimensional regularization 
scheme within n = 4 — 2e dimensions. For the leptogluon mass and field strength renormalization, these quantities 
read: 



5Z U 



47T 



C A [5Je B Q (ml , mj, , 0) - 4 m? 8 3£e B' (mf g , mf g , 0) - 2] 



5m l8 = [2-3ReB (m£,mf 8 ,0)]. 



(24) 
(25) 



Unsurprisingly, the above expressions are identical to the gluon-mediated correction to the gluino field and mass 
renormalization constants. For the gluino case, however, we rely as well on the additional SUSY-QCD contributions 
triggered by the squark interchange. 

Finally, the explicit analytical expressions for the different UV counter terms involving the leptogluon interactions 
5 £ we quote in Table In] in terms of the field, mass and strong coupling renormalization constants. 
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